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Abstract 

No. It is simply not plausible that cosmic acceleration could arise 
within the context of general relativity from a back-reaction effect of in- 
homogeneities in our universe, without the presence of a cosmological con- 
stant or "dark energy." We point out that our universe appears to be 
described very accurately on all scales by a Newtonianly perturbed FLRW 
metric. (This assertion is entirely consistent with the fact that we com- 
monly encounter 5p/p > 10 30 .) If the universe is accurately described by 
a Newtonianly perturbed FLRW metric, then the back-reaction of inho- 
mogeneities on the dynamics of the universe is negligible. If not, then it 
is the burden of an alternative model to account for the observed prop- 
erties of our universe. We emphasize with concrete examples that it is 
not adequate to attempt to justify a model by merely showing that some 
spatially averaged quantities behave the same way as in FLRW models 
with acceleration. A quantity representing the "scale factor" may "accel- 
erate" without there being any physically observable consequences of this 
acceleration. It also is not adequate to calculate the second-order stress 
energy tensor and show that it has a form similar to that of a cosmological 
constant of the appropriate magnitude. The second-order stress energy 
tensor is gauge dependent, and if it were large, contributions of higher 
perturbative order could not be neglected. We attempt to clear up the 
apparent confusion between the second-order stress energy tensor arising 
in perturbation theory and the "effective stress energy tensor" arising in 
the "shortwave approximation." 
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1 Introduction 



The apparent acceleration of our universe is one of the most striking cosmological 
observations of recent times. In the context of Friedmann-Lemaitre-Robertson- 
Walker (FLRW) models in general relativity, the acceleration of our universe 
would require either the presence of a cosmological constant or a new form of 
matter ("dark energy") with large negative pressure. However, since there does 
not appear to be any natural explanation for the presence of a cosmological 
constant of the necessary size, nor does there appear to be any natural candidate 
for dark energy, it is tempting to look for alternative explanations. In recent years, 
there have been at least two approaches that have attempted to account for the 
observed acceleration of our universe within the framework of general relativity 
as being a consequence of deviations of our universe from exact FLRW symmetry, 
without invoking the presence of a cosmological constant or dark energy. 

One approach notes that the mass density of our universe is, in fact, extremely 
inhomogeneous on scales much smaller than the Hubble radius. In order to get 
an effective homogeneous, isotropic universe, one needs to average and/or smooth 
out the inhomogeneities on some appropriate choice of spatial slicing. In such an 
"averaged" (or effective) FLRW universe, one can then define "effective cosmolog- 
ical parameters" [1, 2, 3]. One then finds that the equations of motion for these 
effective cosmological parameters differ, in general, from the equations satisfied 
by these parameters in FLRW models. If they differ in a way that corresponds 
to adding a cosmological constant (or dark energy) of the right magnitude, then 
one may hope to have explained the acceleration of our universe in the context 
of general relativity, without invoking the presence of a cosmological constant or 
dark energy [4, 5, 6, 7, 8, 9, 10]. 

A second approach has attempted to account for the acceleration of our uni- 
verse as a back-reaction effect of long-wavelength cosmological perturbations. 
Here one constructs an "effective energy-momentum tensor" for these perturba- 
tions from second-order perturbation theory, and adds this as a source term in 
Einstein's equation [11, 12, 13, 14]. If this effective energy-momentum tensor has 
a form similar to that of a cosmological constant (or, at least, provides a negative 
pressure term) and is of the appropriate magnitude, then one may hope that this 
could explain the acceleration of our universe. 

Given that our universe appears to be very accurately described by a FLRW 
model, with very small deviations from homogeneity and isotropy, it would seem 
extremely implausible that cosmologically important effects could result from the 
second (or higher) order corrections produced by these small departures from a 
FLRW model. The main purposes of this paper are to explain this point with 
somewhat more precision and to point out some significant flaws in the arguments 
that have been made in the context of the above two approaches. In particular, 
we emphasize that one cannot justify a model by merely showing that spatially 
averaged quantities behave the same way in FLRW models with acceleration; 
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rather one must show that all of the predictions of the model are compatible 
with observations. We illustrate this point by showing that spatial averaging can 
yield "acceleration" for the case of a universe that consists of two disconnected 
components, each of which is decelerating! We also show that spatial averaging 
can yield acceleration for suitably chosen slices of Minkowski spacetime. With 
regard to the back-reaction effects of long- wavelength perturbations, we note that 
the effective energy-momentum tensor is highly gauge dependent, as has previ- 
ously been pointed out by Unruh [15]. Even in the long- wavelength limit, we 
show that one can get essentially any answer one wishes for the effective energy- 
momentum tensor, even though one cannot get any new physical phenomena 
beyond those already present in FLRW models. We also emphasize the differ- 
ence between the use of an effective energy-momentum tensor for gravitational 
perturbations in the context of second-order perturbation theory and the use 
of a similar effective energy-momentum tensor in the context of the "shortwave 
approximation" [16, 17]. The former is highly gauge dependent (and, thus, not 
easily interpreted) and must be "small" if higher order perturbative corrections 
are to be neglected. The latter is essentially gauge independent and need not 
be "small." However, the "shortwave approximation" clearly is not valid for 
analyzing long-wavelength cosmological perturbations. 

In the next section, we point out that our universe appears to be very accu- 
rately described on all scales by a Newtonianly perturbed FLRW metric, despite 
the presence of large density contrasts. If this is correct, then higher order cor- 
rections to this metric resulting from inhomogeneities would be negligible. Thus, 
any model that attempts to explain the acceleration of our universe as a con- 
sequence of higher order effects of inhomogeneities will have to overcome the 
seemingly impossible burden of explaining why the universe appears to be so well 
described by a model that has only very small departures from a FLRW metric. 1 
In section 3, we illustrate that spatial averaging can produce an entirely spurious 
"acceleration" that is not associated with any physical observations. In section 4, 
we emphasize the distinction between second-order perturbation theory and the 
shortwave approximation, and we analyze the gauge dependence of the effective 
energy-momentum tensor arising in second-order perturbation theory. 

2 The Newtonianly Perturbed FLRW Metric 

By a Newtonianly perturbed FLRW metric, we mean a metric of the form 

d~s 2 = -(1 + 2^)dt 2 + a 2 {t)(l - 2^) lij dx i dx j . (1) 

1 For example, rotation of the cosmic matter may produce acceleration effects [5] but these 
acceleration effects must be negligible in view of the observed isotropy of our universe [18, 19, 
20]. 
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where 7^ denotes the metric of a space of constant curvature (3-sphere, flat, or 
hyperboloid) , and \I/ satisfies 

< ^D^D^ , (D'^D^) 2 < (D i D j m)D i D j m , (2) 

where Di denotes the derivative operator associated with 7^. Suppose that the 
stress-energy content of this spacetime consists of fluid components that are very 
nearly homogeneously and isotropically distributed but may have arbitrary equa- 
tion of state (such as radiation, dark energy, and/or a cosmological constant) to- 
gether with components that may be very inhomogeneously distributed on scales 
small compared with the Hubble radius but are nearly pressureless and move 
with velocity much smaller than light relative to the Hubble flow (such as or- 
dinary matter and dark matter). The stress-energy of the smoothly distributed 
components take the form 

T i S) ~ P {s \t)dt 2 + P {s \t)a 2 (t)- fij dx l dx j , (3) 

where P^ = P^ s \p^) is arbitrary, whereas the inhomogeneously distributed 
components have a stress-energy of the "dust" form 

T^W m) (M^ 2 - (4) 

If one plugs the metric form eq. (1) into Einstein's equation and uses eqs. (2), 
the spatial average yields the usual FLRW equations for the scale factor a with 
stress-energy source consisting of the sum of eq. (3) and the spatial average of 
eq. (4), 

3 (^) 2 = W (5) 

3- = -^( p W+pM + 3PW) , (6) 
a 2 

where p denotes the spatial average of p (taken on a t = const, time slice with re- 
spect to the underlying FLRW metric with K — ±1, 0). The dominant remaining 
terms in Einstein's equation then yield 

lA (3) * = y5 P , (7) 

where A/ 3 ) = Y^DiDj, and where 5p = p^ — p^ denotes the deviation of the 
density from the spatial average. In the following we assume the metric 7^ to be 
a flat spatial metric, i.e., K = 0. 

It should be emphasized that the above discussion does not constitute a deriva- 
tion that a metric form, eq. (1), is a good approximation to a solution to Einstein's 



1*1 < 1, 



~dt 
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equation for stress-energy eq. (3) and (4) when eqs (5)-(7) hold. Rather, all that 
has been shown is that if one assumes a metric of the form eq. (1) with eqs. (2) 
holding and with the stress-energy given by eqs. (3) and (4), then eqs. (5), (6) 
and (7) must hold. The above discussion is thus analogous to the usual textbook 
"derivations" of the ordinary Newtonian limit of general relativity, where one also 
postulates a spacetime metric of a suitable "Newtonian form" and assumes that 
the matter distribution is approximately of the form (4) (see, e.g., section 4.4a of 
[21]). One can derive the ordinary Newtonian limit more systematically by con- 
sidering one-parameter families of solutions to Einstein's equation with suitable 
limiting properties (see, e.g., [22] and [23] and references cited therein). It would 
be more difficult to provide an analogous analysis here, but we see no reason to 
doubt that eq. (1) is a good approximation to a solution to Einstein's equation 
when eqs (5)-(7) hold, provided, of course, that conditions (2) are satisfied. 

We now assert that the metric, eq. (1), appears to very accurately describe 
our universe on all scales, except in the immediate vicinity of black holes and 
neutron stars. The basis for this assertion is simply that the FLRW metric 
appears to provide a very accurate description of all phenomena observed on large 
scales, whereas Newtonian gravity appears to provide an accurate description of 
all phenomena observed on small scales. The metric (1) together with eqs (2), 
predicts that large scale phenomena will be accurately described by a FLRW 
model, whereas the metric (1) together with eq. (7) predicts that Newtonian 
gravity will hold in regions small compared with the Hubble radius in which the 
stress-energy (4) dominates over (3) (see [24]). 

Note that the validity of eq. (1) for accurately describing phenomena on small 
scales holds despite the fact that the density contrast of matter is commonly quite 
large [25, 26] 

(8) 

P 

Indeed, for the solar system, galaxies, clusters of galaxies, we can estimate re- 
spectively, Sp/p ~ 10 30 , ~ 10 5 , ~ 10 2 ^> 1. Nevertheless, in all of these cases, we 
have ^ ps 1CT 6 ~ 1CT 5 <C 1, and the other conditions appearing in eq. (2) also 
hold. Even for neutron stars, \I> ps 1CT 1 , so the metric eq. (1) is probably not too 
bad an approximation even in the vicinity of neutron stars. 

The key point of this section is that if our assertion is correct that a metric 
of the form of eq. (1) accurately describes our universe, and if it also is true that 
conditions (2) hold, then the nonlinear correction terms occurring in eqs. (5) and 
(6) are negligibly small. It therefore is manifest that nonlinear corrections 2 to 

2 We take this opportunity to comment upon one misconception related to the validity of 
the Newtonianly perturbed FLRW metric. It is commonly stated that when Sp/p 3> 1, one 
enters a "nonlinear regime." This might suggest that the validity of eqs. (1) and/or (7) would 
be questionable whenever Sp/p >• 1. However, this is not the case; the proper criteria for the 
validity of the metric eq. (1) are conditions (2), not Sp/p <C 1. It is true that nonlinear effects 
become important for the motion of matter when Sp/ p 3> 1. This follows simply from the fact 
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the dynamics of the universe will be negligible, i.e., there will be no important 
"back-reaction" effects of the inhomogeneities on the observed expansion of the 
universe on large scales. In particular, accelerated expansion cannot occur if 
the smoothly distributed matter satisfies the strong energy condition. However, 
our assertion that the metric, eq. (1), very accurately describes our universe is 
merely an assertion, and we cannot preclude the possibility that other models 
(e.g., with large amplitude, long- wavelength gravitational waves or with matter 
density inhomogeneities of a different type) might also fit observations. Our 
main point of this paper, however, is that if one wishes to propose an alternative 
model, then it is necessary to show that all of the predictions of this model are 
compatible with observations such as the observed redshift-luminosity relation 
for type la supernovae and the various observed properties of the cosmological 
microwave background (CMB) radiation. As we shall illustrate in the next two 
sections, it does not suffice to show merely that the spatially averaged scale factor 
behaves in a desired way or that an effective stress-energy tensor is of a desired 
form. 

3 Cosmic Acceleration via Averaging 

The type of spatial averaging in the context of Newtonianly perturbed FLRW 
models that was done to derive eqs. (5) and (6) above is not problematical. The 
metric very nearly has FLRW symmetry, so there is a natural choice of spatial 
slices on which one can take spatial averages. Since \l/ <C 1, it makes negligible 
difference if one uses the spacetime metric (1) or the corresponding "background" 
FLRW metric (i.e., eq. (1) with ^ set equal to zero) to define the averaging. 

If one has a metric that does not nearly have FLRW symmetry, one can, of 
course, still define spatial averaging procedures. However, these will, in general, 
be highly dependent on the choice of spatial slicing, and the results obtained from 
spatial averaging need not be interpretable in a straightforward manner. We now 
illustrate these comments with concrete examples. 

For simplicity and definiteness, we consider an inhomogeneous universe with 

that self-gravitation is a nonlinear effect, and self-gravitational effects on the motion of matter 
cannot be ignored when Sp/p > 1. But this does not mean that one must include nonlinear 
corrections to the metric form, eq. (1), or to eq. (7) in order to get a good approximation to 
the spacetime metric. Indeed, if one is trying to describe the solar system in the context of the 
ordinary Newtonian limit of general relativity, one must include "nonlinear effects" to obtain 
the correct motion of the planets; they would move on geodesies of the flat metric rather than 
the Newtonianly perturbed metric if not for these "nonlinear effects." However, the corrections 
to the spacetime metric of the solar system arising from nonlinear terms in Einstein's equation 
are entirely negligible. 
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irrotational dust. In the comoving synchronous gauge, 3 the metric takes the form 

ds 2 = -dt 2 + qij (t, x m )dx i dx j . (9) 

Let E denote a hypersurface of constant t, let T> denote a compact region of E 
and let be a scalar field on E. The average, {<ft)v, of over the domain DcS 
may be defined by 

{(j>)v=^r [ (10) 
'V Jv 

where Vp denotes the volume of V and dE is the proper volume element of E. 4 
Define the averaged scale factor, a V) by 

a v = {V v yl\ (11) 

We "time evolve" T> by making it be comoving with the dust, i.e., the (comov- 
ing) coordinates of the boundary of T> remain constant with time. Following 
Buchert [1], one then obtains from Einstein's equation the following equations of 
motion for a®, 

~{p)v + Qv, (12) 
K 2 ( P )v-\{n) v - l -Q v , (13) 
together with 

(a 6 v Q v y + a 4 v (a 2 v (K)T>y = 0. (14) 
Here 1Z is the scalar curvature of E and 

Qv = \{{e 2 )v-{e)l)-{a^) v , (is) 

In deriving these equations, it is important to bear in mind that the averaging 
of the time derivative (<p)v = {d(j)/dt)x> of a locally defined quantity <fi differs 
in general from the time derivative of the averaged quantity (4>)' ?, = d(<f)) v /dt, 
since the volume element and Vp may depend on the time t. Indeed, we have 

(fTy = ( 4>)v + - (e)v(c/>)v , (16) 

3 If the universe is filled with irrotational dust, then the comoving synchronous gauge defines 
a natural choice of slicing, namely the slices orthogonal to the world lines of the dust. However, 
for an inhomogeneous universe, this gauge choice typically will break down on timescales much 
shorter than cosmological timescales, due to formation of caustics. For example, synchronous 
coordinates defined in a neighborhood of the Earth would typically break down on a timescale 
of order the free fall time to the center of the Earth, i.e., ~ 1 hour. 

4 For a different type of averaging procedure than that given by eq. (10) and its application 
to cosmology, see e.g., [27, 28]. 
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where 9 denotes the expansion of the world lines of the dust. 

It is immediately seen from (12) that "averaged acceleration" a v > is 
achieved if 



Buchert [29] has discussed cosmological implications of the condition (17). 

A number of authors (see, e.g., Refs. [6, 9, 10]) have sought to account for 
the observed acceleration of our universe by means of inhomogeneous models that 
satisfy eq. (17). However, our main point of this section is that even if our universe 
(or a suitable spatial domain of our universe) satisfies eq. (17) and thus has 
a v > 0, this does not imply that the model will possess any physically observable 
attributes of an accelerating FLRW model. Indeed Nambu and Tanimoto [9] 
have shown that in a cosmological model in which T> can be written as a union 
of regions each of which is locally homogeneous and isotropic, we have 



where ai denotes the locally defined scale factor in the i-ih patch, and, in this 
case, av = (a? + a\ + •••) 1 / 3 . Consider, now, a model where at time t the 
universe consists of two disconnected(l) dust filled FLRW models, one of which 
is expanding and the other of which is contracting. Both components of the 
universe are, of course, decelerating, i.e., a\ < 0, '6,2 < 0. Nevertheless, it is not 
difficult to see from eq. (18) that a© > can easily be satisfied. For example, if 
we take a = a x = a 2 and cti = — d 2 we obtain 



We thereby obtain a very simple model where the universe is accelerating ac- 
cording to the definition eq. (17), but all observers see only deceleration. This 
graphically illustrates that satisfaction of eq. (17) in a model is far from sufficient 
to account for the physically observed effects of acceleration in our universe. We 
see no reason to believe that the spatially averaged acceleration found, e.g., in 
the models of [9, 10] directly corresponds to any physical effects of acceleration 
such as would be observed in type la supernovae data. The only way to tell if a 
model displays physically observable effects of acceleration is to calculate these 
effects. 

As already mentioned above, the averaging procedure defined by eq. (10) also 
has ambiguities both with regard to the choice of time slicing and the choice of 
domain V. One may also artificially produce an averaged cosmic acceleration 
as a result of a suitably chosen time-slicing. To show this explicitly, we give an 
example of accelerated expansion in Minkowski spacetime. We note first that for 



Qv > y(p)o- 



(17) 




(18) 




(19) 
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a general inhomogeneous universe (i.e., with no assumption concerning the form 
of the stress-energy), the equation of motion for a v can be expressed as 

3^ = -<R) V - 6 (^) 2 + ((G ab + \g ab G\)t a t b ) v , (20) 

where G a b denotes the Einstein tensor. Therefore, for any vacuum spacetime, 
if there is a domain V such that —(TZ)v > & {cLv/ a v) 2 i then V describes an 
accelerated expansion insofar as a© is concerned. 

To construct an accelerating V in Minkowski spacetime, we start with two 
hyperboloidal slices, one of which corresponds to an expanding time-slice in the 
Milne chart covering the future of the origin 

ds 2 = -da 2 + a 2 (d£ 2 + sinh 2 £dn 2 ) , (21) 

and the other of which is a similar hyperboloidal slice that is contracting. The 
idea of the construction is to join these two hyperboloids at some radius, smooth 
out the join region, and choose V so that a v = (see Figure 1). Since the 
hyperboloids have negative scalar curvature, we thereby have 1Z < except 
near the radius where the hyperboloids are joined. However, we can show that 
this construction can be done so that the contribution from the join region can 
be made arbitrarily small. Consequently, we obtain 3dv/av = —(7Z)v > 0. 
Details of this construction are given in the appendix. Since Minkowski spacetime 
does not display any physical effects associated with accelerated expansion, this 
example shows quite graphically that "acceleration" as defined by the above 
averaging procedure can easily arise as a gauge artifact produced by a suitable 
choice of time slicing. 

4 Cosmological Back-reaction in the Long- wavelength 
Limit 

A number of authors have considered the back-reaction effects of perturbations 
of FLRW models, particularly with regard to modes whose wavelength is com- 
parable to or larger than the Hubble radius [11, 12, 13, 14]. The basic strategy 
has been to compute the second order terms in Einstein's equation arising from 
these perturbations, thereby obtaining an "effective stress-energy tensor" for the 
perturbations. If the form of this effective stress-energy tensor corresponds to 
that of a positive cosmological constant of the correct magnitude, then one might 
hope to have provided a mechanism for obtaining the observed acceleration of 
our universe as a back-reaction effect of long-wavelength perturbations, without 
the need to introduce a cosmological constant or dark energy. 5 

5 The motivation in [11, 12] was actually to use back-reaction effects to attempt to cancel the 
presence of a large cosmological constant rather than to use back-reaction to directly produce 
acceleration. 
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Figure 1: An accelerating domain T> (thick line) in Minkowski space- 
time (t, r) with the angular coordinates suppressed. The domain V is 
constructed by cutting portions out of the two hyperboloids (thin dashed 
lines), t = V a 2 + r 2 and t = 2\J a? + r 2 . — \/a 2 + r 2 , and joining them 
at r = vq. The matching can be done in a smooth manner, as explained 
in Appendix. The boundary radius r = rp is chosen so that dp = 0. 



We comment, first, that, even without extensive analysis, there is an intrin- 
sic implausibility to this type of explanation. If one considers perturbations of 
wavelength less than the Hubble radius, it is hard to imagine that the pertur- 
bations could be so small that we do not notice any significant deviations from 
homogeneity and isotropy, yet so large that their second order effects produce 
very significant changes to the dynamics of our universe. On the other hand, if 
one goes to the long- wavelength limit, then the perturbation should correspond 
closely to a perturbation to a spatially homogeneous cosmo logical model. But, 
given the severe constraints on anisotropy arising from CMB observations, the 
perturbation should, in fact, correspond to a perturbation towards another FLRW 
model. Thus, one should not be able to obtain any new phenomena (such as ac- 
celeration without a cosmological constant or dark energy) that are not already 
present in FLRW models. 

At least part of the confusion with regard to the calculation of the back- 
reaction effects of cosmological perturbations appears to stem from the fact that 
the notion of an "effective stress-energy tensor" for perturbations arises in two 
quite different contexts, namely (i) ordinary perturbation theory and (ii) the 
"shortwave approximation." We now explain this distinction. For simplicity, we 
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restrict consideration in the following discussion to the vacuum case; cosmological 
perturbations for the Einstein-scalar-field system will be considered later in this 
section. 

Ordinary perturbation theory (see, e.g., section 7.5 of [21]) arises by con- 
sidering a one-parameter family of metrics g a b{ot) that is jointly analytic in its 
dependence on a and the spacetime point. We refer to = g a b(0) as the 
"background metric." Roughly speaking, as a — > 0, g a b(&) differs from g^ by a 
perturbation that becomes of arbitrarily small amplitude but maintains a fixed 
profile. One expands g a b(oi) in a power series in a about a = 



9ab{a) 



= E h a "'* ■ < 22 > 

The perturbation equations for g^ for the vacuum Einstein equation, 

G ab = 0, (23) 

are then obtained by differentiating the Einstein tensor, G a b(a), of g a b(oi) n times 
with respect to a at a = 0. The zeroth order equation is just Einstein's equation 

for 9ab 

G ab [g^]=0. (24) 

The first order equation is 

G^[9 (1) ]=0, (25) 

where denotes the linearized Einstein tensor off of the background metric 
gjjj} . The second order equation is 

G%^] = ~G^], (26) 

where G^[g^} denotes the second-order Einstein tensor constructed from g^ . 

As can be seen from eq. (26), minus the second-order Einstein tensor (divided 
with the gravitational constant), — K~ 2 G^[g^], plays the role of an "effective 
stress-energy tensor" associated with the perturbation g^ in the sense that it acts 

(2) 

as a source term for the second-order metric perturbation g ah . However, this does 
not mean that one can treat —k^ 2 G^ [g^] as though it were a new form of matter 
stress-energy that can be inserted into the right side of the exact Einstein equation 
(23) as opposed to the right side of eq. (26). For one thing, the second-order 
Einstein tensor is highly gauge dependent (as we shall illustrate explicitly below), 
so it is not straightforward to interpret its meaning. The key point, however, is 
that eq. (26) arises only in the context of perturbation theory. If G^[g^} is very 
small, then its effects on the spacetime metric can be reliably calculated from 
eq. (26). But if G^[g^\ is large enough to produce cosmologically interesting 
effects (such as acceleration), then the third and higher order contributions to 
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gabipt) will also be large, and one cannot reliably compute back-reaction effects 
from second-order perturbation theory. 

This situation occurring in perturbation theory contrasts sharply with the 
situation that arises when one uses the "shortwave approximation" [17, 30, 31]. 
Here, one wishes to develop a formalism in which the self-gravitating effects of 
gravitational radiation — and the consequent effects on the spacetime metric on 
scales much larger than the wavelength of the radiation — can be reliably obtained, 
even when these effects are "large." Again, one considers a one-parameter family 
of metrics g a b{(3) that has a continuous limit to the metric g^ = g a b{0). Thus, 
as in ordinary perturbation theory, as (3 — > 0, g a b{(3) differs from by a per- 
turbation of arbitrarily small amplitude. However, one now requires g a b((3) to be 
such that, roughly speaking, as (3 — > 0, the ratio of the amplitude to the wave- 
length of the perturbation goes to a finite, non-zero limit; see [31] for a precise 
statement of what is required in this limit. Thus, in this scheme, the dominant 
terms in Einstein's equation as (3 — > are actually the linear terms in the second 
derivatives of the first order perturbation, which diverge as 1/(3. One thereby 
obtains 

GiM^O. (27) 

The quadratic terms in the first order perturbation are of zeroth order in /3, so 
they make a contribution to the Einstein tensor that is comparable to that of g^ . 
The linear terms in the second-order perturbation also contribute to this order, 
but these contributions can be eliminated by averaging over a spacetime region 
that is large compared with the wavelength of the perturbation. One thereby 
obtains 

G ab [g (0) ] = (-G^[g^]), (28) 

where the brackets on the right side of eq. (28) denote a suitably defined spacetime 
average. It can be shown that (-G { 2[gW\) is gauge invariant in a suitably defined 
sense. We refer to [31] for further details of the derivation and meaning of these 
equations. 

Although eq. (28) is quite similar in form to eq. (26), the meaning and range 
of validity of these equations are quite different. In contrast to eq. (26), it should 
be possible to use eq. (28) to calculate the back-reaction effects of gravitational 
radiation even when these effects are large. The catch, however, is that eq. (28) 
can be used only when the wavelength of the perturbation is much smaller than 
the curvature lengthscale of the background spacetime. Thus, if the universe were 
filled with gravitational radiation of wavelength much smaller than the Hubble 
radius, then it should be possible to use eq. (28) to reliably calculate the back- 
reaction effects of this radiation, even if this radiation is the dominant form of 
"matter" in the universe. However, eq. (28) manifestly cannotbe used to calculate 
the back-reaction effects of long-wavelength perturbations. 

We conclude this section by deriving an explicit formula for the gauge de- 
pendence of the second-order "effective stress-energy tensor" arising in ordinary 
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perturbation theory for long-wavelength scalar-type perturbations of a FLRW 
universe containing a scalar field. By doing so, we will see that one can get 
essentially any answer one wishes for this effective stress-energy tensor by mak- 
ing appropriate gauge transformations. This graphically shows that one cannot 
draw any physical conclusions merely by examining the form of the effective 
stress-energy tensor arising in second-order perturbation theory. 6 
Consider a background flat FLRW universe 

ds 2 = —dt 2 + a 2 {t)^ijdx l dx^ , 

filled with a scalar field whose energy momentum tensor is given by 

T ab = V a 0V 6 - ^a6{V c 0V c + 2C/0)} . (30) 

The unperturbed background equations of motion for a and 0, which are functions 
of only t, are given by 

• dU 
+ 3#0 + — = , 

d(p 

where in these equations and hereafter the dot denotes the derivative with respect 
to t. 

We focus on the scalar-type perturbations. The general form of a scalar-type 
metric perturbation is 

ds 2 = -(1 + 2AS)dt 2 - 2aBS i dtdx i + a 2 {(l + 2# L S) 7ij + 2H T %}dx i dx j , (33) 

and the scalar field perturbation is given by 

= + <50§ . (34) 

Here § denotes a plane wave on flat 3-space with wavevector k, and §; and §y 
are the divergence-free vector and transverse-traceless tensor defined by 

= —DiS , Sij = p {r> i D j - § , (35) 

with Di being the derivative operator associated with the 3-space metric 7^ 
as in (2), and k 2 = k ■ k. Here and in the following, perturbation variables 
are understood as corresponding Fourier expansion coefficients — hence functions 
merely of t — and we omit the index k unless otherwise stated. 

6 Discussion of the characterization of the back-reaction effects of perturbations in terms of 
physical variables can be found in [32, 33]. 
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(29) 



(31) 
(32) 



A - 


-> A - 


-T, 


-T, 

a 


B - 


-> 


- aL + 


H L - 


-> H L 


-HT, 


Hj 1 — 




+ kL, 




5(f) - 


-> 50 


-0T. 





Under infinitesimal gauge transformations of the scalar-type; 

t -> t + T§ , x* -> x* + LS* , (36) 
the perturbation variables, A, B, Hl, Ht, S(f>, change as 

(37) 
(38) 

(39) 

(40) 
(41) 

In particular, it follows from eqs. (38) and (40) that the following combinations 

Xt =KI* t - b )> Xl =-\ Ht > (42) 

change as 

X T ->X T -T, X L ^X L -L. (43) 

Hence, by inspection of eqs. (37), (39) and (41), one can immediately obtain 
gauge-invariant perturbation variables [34, 35]; 

= A-X T , $ = Hl — -X L - HX T , A<\> = 8<j> - <})X T . (44) 

3 

Any scalar-type gauge-invariant perturbation quantity can be expressed as a lin- 
ear combinations of the gauge- invariant variables $, and A(f>, and their time 
derivatives. 

It follows from the linearized Einstein equations that the gauge-invariant vari- 
ables defined above satisfy the following equations [35] 

k 2 r i 

— + = 0, k 2($ - H^) + K 2 (pA(j) =0. (45) 

which correspond, respectively, to the trace-free part of the space-space compo- 
nent and the time-space component of the linearized Einstein equations. For 
k 2 7^ 0, we obtain from eq. (45) the following relations between <3>, \1>, and A<f>; 

$ = ^ = -i + . (46) 

It then also follows from Einstein equations that \& is governed by 



*+(^-|]*+(2|-|]ot + ^ = 0. (47) 
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We therefore have found that all of the scalar-type perturbations variables are 
given in terms of the variables X T) and X L by 



A = ^ + X T , (48) 

B = -aX L --X T , (49) 
a 

H L = -y + HX T + ^X L , (50) 

H T = -kX L , (51) 

5 4- = -^(if + H^+X T . (52) 

The variable \& is gauge invariant and satisfies eq. (47). On the other hand, 
the gauge transformation law (43) implies that the functions Xt and Xl are 
completely arbitrary, i.e., they may be chosen to take any values that one wishes. 
Thus, the specification of Xt and Xl in terms of \1/ corresponds to fixing the 
gauge freedom. For example, the choice 

X T = X L = (53) 

corresponds to the Poisson gauge (or the longitudinal gauge), in which the metric, 
eq. (33), takes precisely the form of eq. (1). Another example is the choice 

X T = -J V(t')dt' + C 1 (k), (54) 

* - '({(^wr^lw)*™' (55) 

where t* denotes some reference time and C\ and C2 are arbitrary constants. 
This choice corresponds to the synchronous gauge, A = B = 0, in which C\ and 
C 2 parameterize the residual gauge freedom in this gauge. 

The second-order effective stress-energy tensor for the Einstein-scalar-field 
system is defined by 

(cff ^--^ ) [^ 1) ]+Ti?[^,a, (56) 

(2) (2) 

where G ab denotes the second order Einstein tensor and T^ b is the similarly 
defined second order contribution to T ab) eq. (30), arising from the first order 
perturbation (5(p, g^). We now calculate the second-order effective stress-energy 
tensor in order to explicitly demonstrate its gauge dependence. It is very conve- 
nient to express ( eS ^T a b in terms of the variables X T) and X L , since any depen- 
dence of ( efC ^T a b on and Xl will explicitly show its gauge dependence. Clearly, 
since ^ eS ^T a b is quadratic in the first order perturbation, it must consist of a part, 
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(eff:*)p a6j faat is quadratic in \I>, a part, ( cS:X ^r a b, that is quadratic in (Xr, -Xi,), 
and a part, ( cff: *' x ^r a fe, containing the "cross-terms" between \1/ and Xt,Xl. The 
quantity <eff:<I '- ) T' a b is gauge invariant, 7 but both ^ efT:X ^T a 6 and ( efT: *' x ^T a ;, are gauge 
dependent. Thus, ^ cS "T a b will be gauge invariant if and only if these latter pieces 
vanish. 

It is easy to verify that both ^ cS:X ^T ab and ^ cS:q, ' x ^T ab are nonvanishing. To see 
explicitly that (c&xy T ab is nonvanishing, it suffices to consider the case where we 
impose the additional restriction 



X T 



-k I* 

Ju a 



x T (t') 



(*0 



dt' , 



(57) 



which ensures that B = 0, thereby considerably simplifying the calculation. We 
also focus attention on the long- wavelength limit. A brute force calculation then 
yields 



k 2 ^T 00 



K 2 {eS:X) T 



9ij 



- (H + 3H 2 )X 2 + (H + 6HH + 12H 3 )X T X T 

3H 2 + 12HH 2 + 3HH + H \ Xl 

2HX T X T - 2(H + 2H 2 )XX - (H + H 2 )X^ 



+ 0(A; 2 )(58) 



-(3H + 22HH + 12H 3 )X T X T 

—HH + AH 2 + 12HH 2 + - 
2 2 



+ 0{k 2 ). (59) 



Thus, even in the case of a pure gauge perturbation, \1> = 0, we can obtain a non- 
vanishing effective stress-energy tensor for long-wavelength perturbations. In- 
deed, since Xt is entirely arbitrary, we see that we can get essentially any answer 
one wishes for ^ eS ^T a b- For example, if one wishes to have a pure gauge perturba- 
tion in which ^"Tab takes the form of a cosmo logical constant, one would merely 
have to solve the second-order ordinary differential equation for Xt that results 
when one equates the right side of eq. (58) to minus the right side of eq. (59). 
This manifestly demonstrates that one cannot derive any physical consequences 
by merely examining the form of the second-order effective stress-energy tensor. 

7 In fact, ^'^^ab is precisely the "effective energy-momentum tensor for cosmological pertur- 
bations" of [11, 12]. However, contrary to the claims of [11, 12], this effective energy-momentum 
tensor is gauge-invariant only in the trivial sense that any gauge-dependent quantity can be 
viewed as gauge invariant once a gauge has been completely fixed. In the variations taken 
in [11, 12] to obtain their effective energy-momentum tensor, Xt and Xl were implicitly as- 
sumed to be independent of 'J, corresponding to the choice of the Poisson (longitudinal) gauge 
Xt = Xl = 0. However, different specifications of Xt and Xl in terms of ^ — i.e., different 
choices of gauge — would lead to different expressions for the effective energy-momentum tensor 
in terms of 5*. 
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Finally, we comment that we derived the above "long-wavelength limit" form 
of ( eS ^T a b f° r scalar-type perturbations by considering perturbations with k 2 ^ 
and then taking the limit as k 2 — > 0. Alternatively, we could have directly 
considered scalar-type perturbations with k 2 = 0. It is not immediately obvious 
that this would give equivalent results, since when k 2 = 0, the quantities §j and 
Sjj do not exist, so B and H T are not defined and the x 1 coordinate freedom 
in eq. (36) does not exist. Furthermore, eqs. (45) become trivial, and therefore 
the relation, eq. (46) need not hold. Thus, it is not entirely straightforward 
to make a physical correspondence between perturbations with k 2 = and the 
k 2 — > limit of perturbations with k 2 ^ 0. Nevertheless, such a one-to-one, onto 
correspondence does exist 8 , and can be explicitly achieved by using the gauge 
freedom available when k 2 ^ to set B = = and using the gauge freedom 
available when k 2 = to set A = —H^. Further discussion of the relationship 
between perturbations in the long- wavelength limit and exactly homogeneous 
perturbations can be found in Refs. [36, 37, 38]. 

Since scalar-type perturbations with k 2 = manifestly correspond to per- 
turbations to other FLRW spacetimes, it is clear that one cannot find any new 
physical phenomena that are not already present in FLRW models by studying 
long- wavelength perturbations and dropping all terms that are 0{k 2 ). For ex- 
ample, consider a FLRW model which contains two matter components, such 
as dust and radiation or two scalar fields. In such a model, there exist non- 
trivial gauge-invariant perturbations even in the k — > limit, and implications 
of such perturbations to the cosmological back-reaction problem have been dis- 
cussed in [8, 32, 33]. However, in the k — > limit such perturbations merely 
correspond to perturbations to other FLRW models; in the above examples, they 
would correspond to changing the proportion of dust and radiation or changing 
the initial conditions of the scalar fields. These perturbations cannot give rise to 
any new phenomena — such as physically measurable acceleration — that are not 
already present in exact FLRW models. 

5 Summary 

In this paper, we have argued that the attempts to explain cosmic acceleration 
by effects of inhomogeneities, without invoking a cosmological constant or dark 
energy, are, at best, highly implausible. A Newtonianly perturbed FLRW metric 
appears to describe our universe very accurately on all scales. In this model, 
the back-reaction effects of inhomogeneities on the cosmological dynamics are 
negligible even though the density contrast may be very large on small scales. 

We focused much of our attention on exposing the flaws in two types of at- 
tempts to explain acceleration by effects of inhomogeneities. (i) Starting from 

8 It also is worth pointing out that, although the gauge freedom is different, the effective 
stress-energy tensor for k 2 = perturbations remains highly gauge dependent. 
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an inhomogeneous model, one can obtain an effective FLRW universe by spatial 
averaging. This effective FLRW universe may display acceleration. However, we 
showed explicitly via concrete examples that acceleration of the effective FLRW 
universe may occur in situations where no physically observable effects of ac- 
celeration actually occur, (ii) The back-reaction effects of a perturbation of a 
FLRW universe are described at second order by an effective stress-energy tensor 
constructed from the first order perturbation. In particular cases, this effective 
stress-energy tensor may take the form of a cosmological constant, thereby sug- 
gesting that it could produce acceleration. However, we pointed out that (unlike 
the effective stress-energy tensor arising in the shortwave approximation), the ef- 
fective stress-energy tensor arising in second order perturbation theory is highly 
gauge dependent and must be small in order to justify neglecting higher order 
corrections. We explicitly evaluated the second-order effective stress-energy ten- 
sor for pure gauge scalar-type perturbations of an Einstein-scalar field model, 
and showed that it can take essentially any form one wishes, including the form 
of a cosmological constant. 
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Appendix 

Here we provide some details of the construction of the accelerating domain T> 
in Minkowski spacetime that was described below eq. (20). Let t and r be, 
respectively, the standard time and radial coordinates in Minkowski spacetime. 
Let a > and r > e > 0. Let / be a smooth, monotone decreasing function of 
one variable such that f(x) = 1 for all a; < 1/2, f(x) = for all x > 1. Define 

Hr) = f + l) • (60) 

Then ifj(r) = 1 whenever r < r — e and ip(r) = whenever r > r . Furthermore, 
there exists a constant C > 0, independent of e, such that e\ip'\ < C for all 
< e < r , where ip' = dip /dr. Define 

F(r) = {^(r) - V(-r + 2r )} (^T^ - yVg + a^j + ^r 2 + a 2 . (61) 

Then F is smooth and the hypersurface S defined by t = F(r) also is smooth. For 
r < r — e, S is the expanding hyperboloid t = \/r 2 + a 2 , whereas for r > r + e, S 
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is the contracting hyperboloid t = — Vr 2 + a 2 + 2^/r 2 + a 2 . The local expansion 
rate H = a/ a of X smoothly changes from 1/a to — 1/a in the junction interval 
(ro — e, ro + e). Nowhere does S display an accelerated expansion locally. 

Now let us take our domain T> to be a ball of radius rp on S, where r£> is chosen 
so that dx> vanishes. It is always possible to find such an rx> since, by construction, 
Vx>, hence dx>, is a smooth function of r, and the local expansion rate H = a/ a is 
positive when r < r — e, whereas it is negative when r + e < r. Since the third 
term of eq. (20) vanishes for Minkowski spacetime, if (TZ)v is negative in D, then 
eq. (20) shows an acceleration 3dx>/ax> = — (7£)z> > 0. However, apart from the 
junction region (r — e, r + e), £ is intrinsically a hyperbolic space with a negative 
scalar curvature 1Z = —Q/a 2 . Furthermore, the following calculation shows that 
the contribution to (TZ)v from the junction region can be made negligibly small. 
The induced metric on E is 

ds 2 = (l-F' 2 )dr 2 + r 2 dn 2 , (62) 

so the scalar curvature of the junction region is given by 

r 2 {1 — (F') 2 } V ; 



Using the formula, 



F , = --f£=(1, + er) + 0(e), (64) 
V r o + a 



which is obtained from the properties of ip in the junction interval, one finds 



ro+e rro+e f o _ ( p/\2 a 

drr 2 ^! - (F') 2 TZ = / dr \ 2 1 J = - 4-g- 



r 



= 0(e), (65) 

which can be made arbitrarily small by taking e — > 0. Thus, the contribu- 
tion to (7£)x> from the junction region can indeed be made arbitrarily small, so 
that -(K) v ~ 6/a 2 > 0. Thus, one obtains an accelerating domain D C S in 
Minkowski spacetime through the volume averaging process. 
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